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Abstract: The determination of scalar lepton and gaugino masses is an important 

part of the programme of spectroscopic studies of Supersymmetry at a high energy e^e~ 

ly-v , hnear colhder. In this article we present results of a study of the processes: e^e~ — )• 

^ ! ^R^~R^ e+e" Xi Xi, e"^e" -^ jjf^ /2^ -^ /i+/x~ Xi Xi, e+e~ -^ e^ e^ -^ e+ e~ X2 X2 

^vq \ and e^e~ -^ Df, I'e ^ e^ e~ xt Xi in two Supersymmetric benchmark scenarios at ^/s 



,— b ' = 3 TeV and 1.4 TeV at CLIC. We characterize the detector performance, lepton energy 

C3 ■ resolution and boson mass resolution. We report the accuracy of the production cross 

section measurements and the Br, JIr, Ve, Xi 1 ^-^d Xi mass determination, estimate the 
systematic errors affecting the mass measurement and discuss the requirements on the 
detector time stamping capability and beam polarization. The analysis accounts for the 
CLIC beam energy spectrum and the dominant beam-induced background. The detector 
performances are incorporated by full simulation and reconstruction of the events within 
the framework of the CLICJLD.CDR detector concept. 
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1 Introduction 

One of the main objectives of linear collider experiments is the precision spectroscopy of 
new particles predicted in theories of physics beyond the Standard Model (SM), such as 
Supersymmetry (SUSY). In this article, we study the production of the supersymmetric 
partners of the muon, electron and neutrino in two specific SUSY benchmark points, where 
we assume R~parity conservation within the so-called constrained Minimal Supersymmetric 
extension of the SM (cMSSM). In this model the neutralino (x?) is the lightest supersym- 
metric particle. Table 1 shows the masses and the branching ratios of the supersymmetric 
particles for the two benchmark points PI and P2. 

For both benchmark points the Higgs boson mass is 120 GeV. Smuons are produced in 
pairs through s-channel 7/Z exchange, selectrons and sneutrinos are pair produced through 
s-channel 7/Z exchange or t-channel x? and Xi exchange respectively, see Figure 1. 

The cross sections, the decay modes, and the cross sections times the branching ratio 
of the signal processes are given in Table 2. In the processes e~^e~ — ?■ ^^^^ each i^ decays 
into an ordinary lepton and a XiJ the x? is stable and escapes detection due to its weakly 



Table 1: Benchmark parameters of the considered SUSY model. 
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Figure 1: Slepton production diagrams: scalar muons (a), scalar electrons (b), and scalar 
neutrinos (c). 



interacting nature. Therefore, the experimental signature of these processes is two oppo- 
sitely charged leptons plus missing energy. For the processes e~^e~ 



e; e 



L^L 



e^ e X2 X2 



and e'^e~ — )■ D^De — )• e~^e~XiXi the signature is an e"'"e~pair, four jets, and missing energy. 
Measuring the lepton energy distributions of these four processes allows the determination 
of their production cross sections and of the eji, flf;, z>e, Xi i aiid Xi masses. The aim of 
this study is to: 

• Characterize the detector performance, namely lepton energy resolution, and boson 
mass resolution. 

• Assess the statistical accuracy of the cross section measurements and the mass de- 
termination. 

• Estimate the systematic errors, affecting the mass measurements, related to the event 
selection and the luminosity spectrum knowledge. 

• Set the requirements for the detector time stamping capability and beam polarization. 

The results presented in this article improve and supersede the previous results [1] obtained 
at 3 TeV only. 



Table 2: Signal processes, decay modes, cross sections, and cross sections times branching 
ratio (cr x Br) at ^/s = 3 TeV and 1.4 TeV. 
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Table 3: Background processes, decay modes and cross sections times branching ratio, 
a X Br, without and with preselection cuts, at 3 TeV and 1.4 TeV. 
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2 Event Simulation 



SUSY signal events and SM background events are generated using the WHIZARD pro- 
gram [2], assuming zero polarisation of the electron and positron beams. WHIZARD is 
interfaced to Pythia 6.4 [3] for fragmentation and hadronization. For the generation 
of processes involving super symmetric particles, the SUSY parameters are entered into 
WHIZARD using the Les Houches format [4] . The physics backgrounds simulated for this 
study are listed in Table 3. 

Beamstrahlung effects on the luminosity spectrum are included using results of the 
CLIC beam simulation for the CDR accelerator parameters [5]. There are three sources of 
the centre-of-mass energy spread: the momentum spread in the linac, the beamstrahlung 
which creates a long tail, and initial state radiation (ISR). The first two are collectively 



refererred to as "luminosity spectrum". The luminosity spectrum is obtained from the 
GuineaPig [6] beam simulation; it is used as input to WHIZARD in which initial state 
radiation and final state radiation (FSR) are enabled. Figure 2 shows the y^ distributions 
for the processes e+e~ — )■ /uj/i^, at -y/s = 3 TeV and e'^e~ — )• eje^ , at -y/s = 1.4 TeV. 
Integrated luminosities of 2000 fb~^ and 1500 fb~^ are assumed at 3.0 and 1.4 TeV respec- 
tively. At ^/s = 3 TeV an integrated luminosity of 2000 fb~^ corresponds to ~4 years (1 
year = 10^ s) of run at the nominal CLIC luminosity of 5.9x10^^ cm~^s~^. At 1.4 TeV 
the nominal luminosity is 3.2x10'^^ cm~^s~^. 

The physics background cross sections of the e+e" — )• i~^(^]i processes are very large, see 
Table 3. Taking into account the luminosity assumptions, the simulation and reconstruction 
of the background events would require very large computing and storage resources. To 
optimize the use of these resources preselection cuts are applied after generation of the 
background events. The preselection requires two opposite charged leptons (LI and L2) 
and the following conditions: 

• pt{L1 and L2) > 4 GeV and 10° < 6(11 and L2) < 170° 

• 4° < A(/)(L1,L2) < 176°, pt{L1,L2) > 10 GeV and M{L1,L2) > 100 GeV 

where pT is the transverse momentum, 6 the polar angle of the lepton, A(j){Ll,L2) the 
acoplanarity of the leptons, pt(L1)L2) the vector sum of the px of the two leptons, and 
M(L1,L2) the invariant mass of the two leptons. Table 3 shows the decay modes, and 
the cross section times branching ratio values without and with preselection cuts. For 
the signal samples, these cuts are also applied after full simulation and reconstruction. 
The simulation is performed using the GEANT4-based [7] Mokka program [8] with the 
CLICJLD_CDR detector geometry [9], which is based on the ILD detector concept [10] 
being developed for the ILC. 

3 Event Reconstruction 

Events are subsequently reconstructed using the Marlin reconstruction program [11]. The 
tracking systems of the CLIC detectors are designed to provide excellent momentum mea- 
surement for charged particle tracks. The track momenta and calorimeter data are input 
to the PandoraPFA algorithm [12, 14] which performs particle flow (PFO) reconstruc- 
tion, including particle identification and returns the best estimate for the momentum and 
energy of the particles. 

3.1 Tvi^o Lepton final states 

The energy of the lepton is reconstructed from the momentum of the charged particle track 
and corrected for final state radiation and bremsstrahlung. The energy of photons and e~^e~ 
pairs from conversions within a cone of 20° around the reconstructed lepton direction is 
added to the energy from the track. Figure 3 shows, for the process e+e" — )• e^e^, the true 
and reconstructed lepton energy distributions without (a) and with (b) photon radiation 
correction. For the process e'^e~ — )• p'^fj-^ the photon radiation corrections are much 
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Figure 2: Centre-of-mass energy spectrum for the processes: e+e — )■ /i/^/i^, \/s = 3 TeV 
(a), and e+e~ — )• e^e^ , \/s = 1.4 TeV (b). 

Table 4: Reconstruction efficiency, e/j, without and with 77 — )• hadrons overlaid for the 
different signal processes, at y^ = 3 TeV and 1.4 TeV. The statistical error on these 
efficiencies is ~ 0.5%. 
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smaller. For both processes there is a good agreement between the true and reconstructed 
lepton energy distributions when photon radiation corrections are applied. 

Table 4 shows the reconstruction efficiencies, ep,, for the signal processes. For the 
process e+e" — )• i~^^~^, er is the number of good reconstructed lepton pairs divided by the 
number of generated lepton pairs. A lepton is considered as good when the reconstructed 
lepton matches the generated particle in space within 2°. For the process e"'"e~ — )• iJ-^fJ-]^, 
at 3 TeV and 1.4 TeV, there is an inefficiency of about 2.5%; 2.0% is due to the cut 
on the lepton angle and 0.5% is coming from muon misidentification. For the process 
e+e~ — )■ e^e]^, at 3 TeV, there is an inefficiency of 6.5%; 4.0% is due to the cut on the 
lepton angle and 2.5% is coming from electron reconstruction or misidentification. At 
1.4 TeV the inefficiency is 5.5%; 3.0% is due to the cut on the lepton angle and 2.5% is 
coming from electron reconstruction or misidentification. 

The energy resolution is characterized using: AE / E^^^^, where AE = £^True — -E'rcco, 
£^True is the lepton energy at generator level before final state radiation or bremsstrahlung, 
and £'Reco is the reconstructed lepton energy with photon radiation corrections. 



Figure 5 (a) and (b) show the lepton energy resolution, for the two lepton final state 
processes at ^/s = 3 TeV and without beam induced background, 77 — )• hadrons. The 
resolution is parametrised using the sum of two Gaussian functions Gl and G2; Gl for 
the peak and G2 for the tails. For the muons the r.m.s. of Gl is 1.5 • 10~^ GeV" , and 
the r.m.s. of G2 is 4.9 • 10~^ GeV~ . Only 4.1% of the events are outside of the central 
region; the central region of the distribution is defined within the interval IS.E/E'^^^^ = 
±0.5 • 10~^ GeV~ . The electron energy resolution is described by the Gaussian Gl with a 
very similar r.m.s. as that for muons, 1.4- 10~^ GeV~ , however, even with bremsstrahlung 
recovery, about 30% of the events are outside the central region. These are due to cases 
where final state radiation and bremsstrahlung are not sufficiently well accounted for; the 
tails are reasonably well described by the Gaussian G2 with r.m.s. = 7.7 • 10~^ GeV^ . 

3.2 Tvi^o leptons and four jets final states 

For the processes e^e~ — )• e^e^ and e^e~ — )■ Ugi^e, the parton topology signature required 
is two leptons and four quarks. After the reconstruction of all the particles in the event, 
the jet finder program FastJet [16] is used to reconstruct jets. The jet algorithm used 
is the inclusive anti-kt method [17]; The choice of cylindrical coordinates is optimal since 
the 77 —7- hadrons events are forward boosted, similarly to the underlying events in pp 
collisions for which the anti-kt clustering has been optimised. The R parameter cut value 
is 1 and the minimum jet energy required is 20 GeV at -y/i = 3 TeV and 10 GeV at 1.4 TeV. 
An event is retained if six jets are found and if two of the jets are identified as isolated 
leptons. Table 4 shows the reconstruction efficiencies of both processes, er is the number of 
reconstructed six jet events, with two leptons, divided by the number of generated events 
with two leptons and four quarks. 

Figure 4 shows the electron energy distribution for the processes e^e~ — )• VeT^e (a) and 
e+e~ — )■ e^e^ (b). There is good agreement between the true and reconstructed electron 
energy distributions when photon radiation corrections are applied. 

For the processes with two electrons and four jets, see Figure 5 (c) and (d), despite the 
presence of four jets, the electron energy resolution is consistent with the energy resolution 
obtained for the isolated electrons process, see Figure 5 (b). 

3.3 Reconstruction ^vith beam-induced background 

The creation of electron-positron pairs and the production of hadrons in 77 interactions are 
expected to be the dominating source of background events originating from the interaction 
region [15]. The beam-beam interaction leading to the production of these background 
particles was simulated with the GUINEAPIG program [6]. The average number of 77 
interactions for each bunch crossing is 3.2 at 3 TeV and 1.3 at 1.4 TeV. At 3 TeV the 
pile-up of this background over the entire 156 ns bunch-train deposits 19 TeV of energy in 
the calorimeters, of which approximately 90% occurs in the endcap and 10% in the barrel 
regions. On average, there is 1.2 TeV of reconstructed energy from 77 — t- hadrons that 
are in the same readout window as the physics event. To reduce this energy deposit, px 
and additional timing cuts are applied. The presence of the 77 — t- hadron background sets 
strong requirements for the design of the CLIC detector and its readout. 
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Figure 3: Process e+e — )• e^e^ij at ^/s = 3 TeV: true and reconstructed electron energy 
distributions, without (a), and with photon energy correction (b). 
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Figure 4: True and reconstructed electron energy distribution, with photon energy cor- 
rection, for processes e+e" — )• i''^i'~ (a), and e+e~ — )■ e^e^ (b) at y/s = 3 TeV. 



To investigate the effect of beam-induced background, the reconstruction software is 
run overlaying particles produced by 77 — )■ hadrons interactions [18]. The 77 — )• hadrons 
event sample was generated with Pythia and simulated. From this sample we randomly 
select for each physics event the equivalent of 60 bunch crossings, assuming 3.2 events per 
bunch crossing at 3 TeV [15] and 1.3 events per bunch crossing at 1.4 TeV. 

The detector hits from these events are merged with those from the physics event 
before the reconstruction. A time window of 10 nsec on the detector integration time is 
applied for all detectors, except for the HCAL barrel for which the window is 100 nsec. 
After particle reconstruction timing cuts in the range of 1 to 3 nsec are applied in order 
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Figure 5: Lepton energy resolution, at ^/s = 3 TeV, without 77 — )• hadron background, for 

the processes: e+e~ — )• fj^'^fj''^ (a), e+e" — >• e^e]^ (b), e+e~ — >• t'gt'e (c), and e+e~ 

(d). 
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to reduce the number of particles coming from 77 — )• hadrons interactions and to optimize 
the energy resolution. The cut values vary according to the particle type (photon, neutral 
hadron, charged particle), the detector region, (central, forward) and the pr of the particle. 
Table 5 shows the cut values for the tight particle flow (PFO) selection. 

Figures 6 (a) and (b) show the lepton energy resolution, without and with 77 — t- hadron 
background, for the processes e+e~ — )• /i^/i^ and e+e~ — )■ e^e^ respectively; only a 
cut requiring p^ > 4 GeV is applied. The lepton energy resolution is preserved; the 
event selection efficiency is reduced by 1.0% for e+e~ — )• fL^fi'^ and 3.3% for e+e~ — >• 
^%,^~R ' ^^^ Table 4. At 1.4 TeV the 77 — )• hadron background is a factor two lower, no 



Table 5: Tight particle flow (PFO) selection cuts for the 77 — )• hadron background sup- 
pression. The cut value depends on the particle id, on the detector region and on the px- 
The same cuts are used at ^/s = 3 TeV and ^/s = 1.4 TeV. 
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selection inefficiency is induced for the process e^e~ — )• /i^A^; an inefficiency of 1.5% is 
induced for the process e^e~ — )• e^e^. In final states with four jets and two leptons, the 
background from 77 — )-hadrons cannot be removed using a similar pT cut, as this would 
significantly degrade the jet energy reconstruction. Figure 6 (c) shows the bias in the 
reconstructed electron energy when the 77 — )-hadron background is included. This bias is 
due to additional background particles being associated with the electron in the attempt 
to account for FSR and bremsstrahlung. Without PFO cuts, the energy resolution is not 
degraded but the central value is shifted. Figure 6 (d) shows the lepton energy resolutions 
without and with 77 — )■ hadrons overlaid after tight PFO selection cuts. The cuts restore 
the central value and preserve the energy resolution, but reduce the reconstruction efficiency 
eR by 6%, see Table 4. 

Figure 7 shows, for the process e^e~ — )• VeT^e at 3 TeV, the W boson mass distribution 
without and with overlaid background: without PFO cuts (a) and with tight PFO selection 
cuts (b). The tight selection cuts give a similar mass distribution as the one obtained with- 
out overlaid background. To estimate the mass resolution degradation. Figure 8 shows the 
W boson mass distribution fit, for the process e^e~ — )• z^gt'e without overlaid background 
(a) and with overlaid background and tight selection cuts (b). The mass distributions are 
fitted with a Breit-Wigner convoluted with two Gaussians, one Gaussian takes into account 
the resolution in the peak, the second the tails. The most probable mass value is fixed as 
well as the natural width of the W . The width of the peak convoluted Gaussian is 4.1 GeV 
without overlaid background, it increases to 4.7 GeV with overlaid background and tight 
PFO selection cuts. The fraction of events in the peak gaussian is 90% without overlaid 
background and 89% with overlaid background. Figure 9 (a) shows the W and H boson 
mass distributions for the processes e^e~ — t- v^Ve and e^e~ — t- e^e^. The distributions 
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Figure 6: Lepton energy resolution, witiiout and with overlaid background at 
^/s = 3 TeV, for the processes: e+e" — )• /i^A^, pt > ^ GeV; no PFO selection (a), 
e+e~ — )• e^e^, pt > 4 GeV; no PFO selection (b), e^e~ — )• t'et'e; no PFO selection (c), 
and e^e" — ?• i'ei'e^ tight PFO selection (d). 



which correspond to an integrated luminosity of 2000 fb^^ are fitted with two Breit-Wigner 
functions. The mass distribution of the Higgs boson is broader than the W one, due to a 
10% background component from Z boson decays and due to semi-leptonic heavy flavour 
decays in the H ^ bb process. In this analysis, no flavour tagging is applied. Figure 9 
(b) shows the boson mass distributions for all inclusive SUSY processes with four jet final 
states [13]. It illustrates that adding 6-tag information in the analysis would improve the 
separation of W and H final states. 
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Figure 7: W boson mass distribution, at ^/s = 3 TeV, for the process e~^e — )■ D^De without 
and with overlaid background: without PFO cuts (a) and with tight PFO selection (b). 
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Figure 8: W boson mass distribution fit, at ^/s = 3 TeV, for the process e~^e — )• Pgi^e 
without overlaid background (a) and with overlaid background and tight selection cuts (b). 

4 Event Selection 

All signal processes have two undetected x^'s in the final state. Therefore, the main char- 
acteristics of these events are missing energy, missing transverse momentum and acopla- 
narity. Despite this signature, the large Standard Model backgrounds make the analysis 
rather challenging. To distinguish signal events from background events the following set 
of discriminating variables is used: 
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Figure 9: W and H boson mass distributions, at y^ = 3 TeV, for 
the processes: e~^e~ — )■ Pg'^e and e~^e~ — )■ e^e^ — )• 2 leptons and 4 jets (a), and 
e^e~ —7- inclusive SUSY — )• 4 jets (b). 

• dilepton energy E{L1) + E{L2), 

• vector sum pt{L1, L2) of the two leptons, 

• algebraic sum pt{L1) + pt{L2) of the two leptons, 

• dilepton invariant mass M(L1,L2), 

• dilepton velocity /3(L1,L2), 

• cos0(Ll, L2); 0{L1, L2) is the polar angle of the vector sum of the two leptons, 

• dilepton acollinearity it — 02 — 9i, 



dilepton acoplanarity vr 



n, 



• dilepton energy imbalance A = \E{Ll) - E{L2)\/\E{L1) + E{L2)\, 

where LI and L2 are the two leptons. For illustration, Figure 10 (a) shows for the process 
e^e~ — )• ji^jl~^, the normalized distributions of some of the observables, for signal and 
background events, namely the dimuon energy, the vector sum of the pT of the leptons, 
the algebraic sum of the pT of the leptons the dimuon invariant mass, the acolinearity and 
the polar angle of the vector sum of the leptons. 

The event selection proceeds as follows. The signal and background samples are split 
into two equal data size samples called "Monte Carlo" and "Data". The events of each 
sample are weighted such that the samples correspond to the same integrated luminosity. 
Then the Boosted Decision Trees (BDT) method from the multivariate analysis toolkit. 
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Figure 10: Discriminating variables used to separate signal and background events for 
the process e~^e~ — )• /i^/i^ at ^/s = 3 TeV: dimuon energy (a), dimuon pT vector sum (b), 
dimuon px algebraic sum (c), dimuon invariant mass (d), acolinearity (e) and polar angle 
of the vector sum of the leptons (f). 



TMVA[19], is used to implement the event selection. Firstly the discriminating variables 
of the Monte Carlo sample are input to the BDT method which trains the BDT probability 
classifier and computes the weights allowing to distinguish signal from background. Next 
the weigths are used to the evaluate the "Data" sample, computing for each event a prob- 
ability value allowing to rank the events to be signal or background-like. The cut value is 
chosen to optimise the significance Smc/V^mc + Bj^jc versus the signal efficiency and the 
background rejection; Smc and B]\jc are the number of signal and background events of the 
MC sample. The cross section and the masses are determined after background subtraction 
and efficiency correction; the errors on the masses depend on \/Sdata + B^ata + Bmc- A 
stronger BDT cut reduces slightly the significance but decreases significantly the errors on 



the masses. Figure 11 shows for the process e~^e 



f? R 



^ = 1.4 TeV, the stacked 
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Figure 11: Process e^e — )• e]^e^ at ^/s = 1.4 TeV: electron energy distribution for signal 
and background events with loose BDT cut (a), and with optimized BDT cut (b). 



electron energy distribution for signal and background events with a loose BDT cut (a), 
and with an optimized BDT cut (b). At 3 TeV the BDT selection efficiency is 95% for 
the dimuon events, 90% for the dielectron events and 94% for the dielectron and four jet 
events. At 1.4 TeV the efficiency is 90% for the dimuon events, 80% for the dielectron 
events, and 90% for the dielectron and four jet events. 



5 Slepton and Gaugino Mass Determination 

After the final selection, the slepton, neutralino or chargino masses are extracted from the 
position of the kinematic edges of the lepton energy distribution, a technique first proposed 
for squarks [20], then extensively applied to sleptons [21]: 



m 



i± 



4n 



{Eh - ELf 
{Eh + ELf 



1/2 



and m^o or m-± 



m 



l± 



1 



2{Eh + El) 



where El and Eh are the low and high edges of the lepton energy distribution 



1/2 



(5.1) 



Eh, l 



2 
X? 



1± 



,"^L 



m- 



e±, 



(5.2) 



The masses are determined using a three-parameter fit to the background subtracted 
energy distribution, with a^±, m^± and m^o or rn~± as parameters. The background 
subtraction is done using the "Monte Carlo" event sample used to train the classifier. The 
fit is performed with the Minuit minimization package [22]. The fit function is: 

f{E)= I LEff{Vs)- I _ U{ai±,m'^±,m^Q,y/s,E -t) ■ D{t) d^/s dr (5.3) 



Y -^r) 



El{V~s) 
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LEff{yfs) is the effective luminosity function, LEff{yfs) = L{yfs) ® ISR{yfs) (g) a"^±(y^). 
L{^/s) is the luminosity spectrum prior to initial state radiation (ISR), ISR{-^fs) is the -y/s 
variation due to ISR and a^± i^/s) is the slepton cross section. [/ is a uniform distribution 
of E, and depends on the process cross section a^±, the slepton and gaugino masses and 
y^ ; the boundaries El, Eh of U are given by 5.2. D is the detector resolution function 
obtained from the fits shown in Figure 5. Figure 12 shows, for the processes e~^e~ — )• 
e^e^ (a) and e^e~ — >■ D^Df, (b) at ^/s = 3 TeV the lepton energy distributions and fit 
results. Table 6 shows the values of the measured slepton cross sections, slepton masses, 
and gaugino masses at ^/s = 3 TeV, assuming 2 ab~^ of integrated luminosity. For the 
process e~^e~ — )■ e^ e^ — )• e^ e~ X2 X2^ the cross section is determined from the fit to the 
boson mass distribution. Figure 9. Table 7 shows the results at 1.4 TeV, assuming 1.5 ab^^ 
of integrated luminosity. 
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Figure 12: Lepton energy spectrum and fit results, for the processes: e+e — )• e^e^ (a) 
and e'^e" — )• t'et'e (c) at -y/i = 3 TeV. 



Table 6: Values of cross sections, slepton and gaugino masses, and statistical accuracies 
assuming an integrated luminosity of 2 ab~^ at ^/s = 3 TeV. 



Process 


Decay Mode 


a 


mi 


m^o or m-± 

Al Xl 






fb 


GeV 


GeV 


e+e^ -^ fl^fi]^ 


,,+ ,,-c,0.<-0 


0.73 ± 0.02 


1011.9 ± 4.9 


342.7 ± 9.7 


e+e- ^ e+e^ 


e e XiAi 


6.23 ± 0.04 


1012.8 ± 2.9 


345.7 ± 4.6 


e+e- -^ elel 


e+e-x^X^ 


2.77 ± 0.20 






e+e- —7- PgPe 


e+e-x^X? 


13.27 ± 0.23 


1094.0 ± 2.3 


644.0 ± 3.6 
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Table 7: Values of cross sections, slepton and gaugino masses, and statistical accuracies 
assuming an integrated luminosity of 1.5 ab~^ at ^/s = 1.4 TeV. 



Process 


Decay Mode 


a 


m^ 


m^o or m-± 






fb 


GeV 


GeV 


e+e" -^ jx^fl~^ 


/i+^-X?X? 


1.51 ± 0.03 


559.1 ± 0.4 


357.1 ± 0.7 


e+e- ^ e+e;^ 


e+e-x?X? 


5.99 ± 0.05 


557.9 ± 0.6 


356.1 ± 0.9 


e+e" — )• i>et'e 


e+e-x^X^ 


5.13 ± 0.19 


644.5 ± 2.2 


488.8 ± 1.1 



6 Systematic Uncertainty related to the event selection 

For the event selection described in section 4 the signal sample, used to train the classifier 
allowing to distinguish signal events from background events, was generated with the same 
slepton and gaugino masses as the data sample. With real data the masses are unknown. 
In this section we describe the procedure allowing to determine the masses and assess the 
error on the masses introduced when the MC masses are different from the true masses; 
the evaluation is done for the process e+e" — )■ /iJ/U^ at 1.4 TeV. 

Firstly signal events for lower smuon and neutralino masses are simulated and recon- 
structed; the smuon and neutralino masses are 459 GeV and 257 GeV respectively, that is 
to say, data masses - 100 GeV. These events are used to train a classifier in which three 
variables are removed, namely the dilepton energy, dilepton velocity and dilepton energy 
imbalance. These variables are most correlated with the masses. The 6 variables classifier 
is then used to select the events. Figure 13 (a) shows the the stacked muon energy distri- 
bution for signal and background events selected with the 6 variables classifier trained with 
masses lower by 100 GeV. The energy distribution of the MC training sample is the black 
dotted line; the energy distribution of the selected signal data sample is the black full line; 
and the energy distribution of the signal data sample without selection is the black dashed 
line. 

Next signal events for larger smuon and neutralino masses are simulated and recon- 
structed. The smuon and neutralino masses are 659 GeV and 457 GeV respectively, this is 
to say, data masses + 100 GeV . These events are used to train the classifier which is then 
used to select the events. Figure 13 b) shows the the stacked muon energy distribution for 
signal and background events selected with a classifier trained masses with masses larger 
by 100 GeV. 

The signal energy distribution of the sample trained with larger masses is obviously 
biased, nevertheless the end points are visible and have similar values as the ones of the 
sample trained with lower masses. A rough estimation leads to El=80 GeV and -E'//=340 
GeV; with these values and y(s)=1.4 TeV formula 5.1 leads to a fl mass of 549.8 Gev 
and a Xi mass of 347.7 GeV, these values are about 10 GeV lower than the true mass 
values. These mass values are then used to simulate and reconstruct a new signal sample 
which is used to train the classifier with all 9 variables. Figure 14 a) shows the the energy 
distribution of the data sample selected with the 9 variables classifier trained with masses 
lower by 10 GeV and the fit result. Figure 14 b) shows the the energy distribution of the 
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Figure 13: Process e^e~ — )• /ijj/ij^ at ^/s =1.4 TeV: stacked inuon energy distribution 
for signal and background events with BDT trained wither lower signal masses (a), and 
with BDT trained wither larger signal masses (b). 

data sample selected with a the 9 variables classifier trained with the true masses and the 
fit result. The smuon and neutralino masses determined after selection of the events with 
the classifier trained with masses lower by 10 GeV are lower by 1.0 GeV and 1.3 GeV 
respectively; the values are statistically compatible 



7 Systematic Uncertainty related to the Luminosity Spectrum 

The beam energy is derived from the beam deflection measurement using high precision 
beam position monitors (BPM) pairs placed before and after the first dipole in the energy 
collimation section. This setup provides a relative energy resolution better than 0.04% 
[23]; therefore the impact on the slepton and gaugino masses is considered as negligible. 

In this section the systematic uncertainty on the slepton and gaugino masses, re- 
lated to uncertainties in the the knowledge of the luminosity spectrum, is investigated. 
The assessment is done at 3.0 TeV where the beamstrahlung is largest. As can be seen 
from equation (5.3), the slepton and gaugino masses depend on the effective luminosity 
function LEff{y/s, p) = L{y/s, p) ^ ISR{^/s) ^ aj;±{y/s). The details about the method 
used to reconstruct the luminosity spectrum L{y/s, p) using Bhabha events are reported 
in [24]. The luminosity spectrum is parametrized with a function F{xi,X2,p) where 
xi^2 = 2-Ei^2/\/s; -£'1,2 is the energy of the e"'"e~ particles before ISR; the vector p has 19 
parameters. The model takes into account the longitudinal boost, the correlation between 
the two particle energies and accounts for asymmetric beams. A fit of F(xi,X2,'p) to 
the Bhabha events using the energy and the acollinearity of the outgoing e'^e~ particles 
allows to determine the parameters p of the luminosity function and their errors. The 
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Figure 14: Process e+e~ — t- /i^/ij!j at \/i = 1.4 TeV: muon energy distribution and fit 
results for events selected with a 9 variables classifier trained with signal masses lower by 
10 GeV (a), for events selected with a 9 variables classifier trained with the true signal 
masses (b). 

parameters were determined at 3 TeV, using 2.2 • 10^ events and taking into account the 
e~^e~ energy resolution. 

To estimate the systematic error on the masses due to the luminosity spectrum, 
the mass fit is performed 38 times. Prior to each fit the effective luminosity spectrum 
LiEff{\fs, ~p) is computed; one parameter pi is changed to pi + ^ or pj — ^ and all other 
parameters are kept to their nominal value. 
The error on the mass from the luminosity is: 



(Tn 






Cij is the correlation matrix obtained from the luminosity spectrum fit and: 
5i = miLEffiVs, f + ^i^)) - miLEffiVs, f - ^i^))- 



(7.1) 



(7.2) 



where m is the the result of the mass fit described in section 5. 

Table 8 shows the values of the slepton and gaugino masses, the corresponding statisti- 
cal uncertainty, and the systematic errors from the knowledge of the shape of the luminosity 
spectrum. For 2 ab~^ of integrated luminosity, the statistical errors are dominant. 

8 Polarization 

Beam polarization is very helpful in the study of SUSY processes both to improve the signal- 
to-background ratio and as an analyzer [25] , in particular to establish the chirality of the 
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GeV 
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(lumi)% 


GeV 


(stat)% 


(lumi)% 


e+e" -^ fl^flJi^ 


1011.5 


0.45 


0.02 
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2.8 


0.06 


G C r CdOd 


1011.8 
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0.02 


341.2 


0.80 


0.06 


e+e~ ^ i>ei^e 


1093.6 


0.19 


0.02 


643.6 


0.54 


0.03 



Table 8: Slepton and gaugino masses, statistical and systematic uncertainties, from the 
knowledge of the shape of the luminosity spectrum (lumi) , at ^/s = 3 TeV. 

sleptons. Table 9 shows the signal cross sections for different electron and positron beam 
polarization conditions. Running with left polarized electron beam would establish the 
chirality of the selectron which decays into two leptons and of the selectron and sneutrinos 
which decay into two leptons and four jets. Running with right polarized electron beam 
would increase the cross sections of the in processes and reduce some of the backgrounds. 

Table 9: Signal processes cross sections {a), for different electron, positron beam polar- 
ization conditions, at \/i = 3 TeV. 



Process 


e+e -^ jl^i^^ 


e^e ^^R^R 


e+e ^e+^L 


e+e —7- i'ei'e 




<m 


^m 


<fh) 


<m 


beam polarization 










e~ 


none , e"*" mone 


0.72 


6.05 


3.06 


13.76 


e~ 


L80%, e+ mone 


0.46 


2.59 


4.78 


21.90 


e~ 


R80%, e+ mone 


0.98 


9.51 


1.34 


5.62 


e~ 


R80%, e+ :L60% 


1.15 


11.40 


1.14 


4.56 



9 Summary 

The accuracy of the slepton and gaugino mass determination and of the process cross 
section measurement in pair produced eji,eL,jJ'R, and u^ processes has been studied at 
CLIC with the CLIC_ILD_CDR detector model for two specific SUSY benchmark scenarios 
at y/s = 3 TeV and 1.4 TeV. The analysis is based on two lepton and two lepton plus four 
jet final states. 

The electron and muon energy resolution and the boson mass resolution are not affected 
by the beam induced background, provided the detectors have timing capabilities of the 
order of 1 nsec allowing for the application of PFO selection cuts. The reconstructed boson 
mass accuracy allows W and light H final states to be distinguished; b tagging improves 
the purity of the W and H samples. 

Slepton cross sections, slepton and gaugino masses can be extracted from the lepton 
energy distributions. At 3.0 TeV, for 2.0 ab~^ of integrated luminosity the relative statisti- 
cal error on the masses is in the range of 0.15 to 0.45% for the sleptons and in the range of 
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0.5 to 2.8% for the gauginos. At 1.4 TeV, for 1.5 ab^^ of integrated luminosity, the relative 
statistical errors, on the slepton and gaugino masses are in the range of 0.1 to 0.2%. 

A major source of smearing of the kinematic edges of the lepton energy spectrum 
is beamstrahlung and ISR. The measurement of the luminosity spectrum with Bhabha 
events, allows a good control of the beamstrahlung. The systematic errors on the slepton 
and gauginos masses due to the knowledge of the luminosity spectrum were estimated. 
At 3.0 TeV for 2.0 ab~^ of integrated luminosity the statistical errors are larger than the 
systematic errors. 
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